Knocking Out Chloroplastic Aldolases/Rubisco Lysine Methyltransferase Enhances Biomass Accumulation in Nannochloropsis oceanica under High-Light Stress.

Rubisco large-subunit methyltransferase (LSMT), a SET-domain protein lysine methyltransferase, catalyzes the formation of trimethyl-lysine in the large subunit of Rubisco or in fructose-1,6-bisphosphate aldolases (FBAs). Rubisco and FBAs are both vital proteins involved in CO2 fixation in chloroplasts; however, the physiological effect of their trimethylation remains unknown. In Nannochloropsis oceanica, a homolog of LSMT (NoLSMT) is found. Phylogenetic analysis indicates that NoLSMT and other algae LSMTs are clustered in a basal position, suggesting that algal species are the origin of LSMT. As NoLSMT lacks the His-Ala/ProTrp triad, it is predicted to have FBAs as its substrate instead of Rubisco. The 18-20% reduced abundance of FBA methylation in NoLSMT-defective mutants further confirms this observation. Moreover, this gene (nolsmt) can be induced by low-CO2 conditions. Intriguingly, NoLSMT-knockout N. oceanica mutants exhibit a 9.7-13.8% increase in dry weight and enhanced growth, which is attributed to the alleviation of photoinhibition under high-light stress. This suggests that the elimination of FBA trimethylation facilitates carbon fixation under high-light stress conditions. These findings have implications in engineering carbon fixation to improve microalgae biomass production.

Mutation mapping of a variegated EMS tomato reveals an FtsH-like protein precursor potentially causing patches of four phenotype classes in the leaves with distinctive internal morphology.

BACKGROUND: Leaf variegation is an intriguing phenomenon observed in many plant species. However, questions remain on its mechanisms causing patterns of different colours. In this study, we describe a tomato plant detected in an M2 population of EMS mutagenised seeds, showing variegated leaves with sectors of dark green (DG), medium green (MG), light green (LG) hues, and white (WH). Cells and tissues of these classes, along with wild-type tomato plants, were studied by light, fluorescence, and transmission electron microscopy. We also measured chlorophyll a/b and carotene and quantified the variegation patterns with a machine-learning image analysis tool. We compared the genomes of pooled plants with wild-type-like and mutant phenotypes in a segregating F2 population to reveal candidate genes responsible for the variegation. RESULTS: A genetic test demonstrated a recessive nuclear mutation caused the variegated phenotype. Cross-sections displayed distinct anatomy of four-leaf phenotypes, suggesting a stepwise mesophyll degradation. DG sectors showed large spongy layers, MG presented intercellular spaces in palisade layers, and LG displayed deformed palisade cells. Electron photomicrographs of those mesophyll cells demonstrated a gradual breakdown of the chloroplasts. Chlorophyll a/b and carotene were proportionally reduced in the sectors with reduced green pigments, whereas white sectors have hardly any of these pigments. The colour segmentation system based on machine-learning image analysis was able to convert leaf variegation patterns into binary images for quantitative measurements. The bulk segregant analysis of pooled wild-type-like and variegated progeny enabled the identification of SNP and InDels via bioinformatic analysis. The mutation mapping bioinformatic pipeline revealed a region with three candidate genes in chromosome 4, of which the FtsH-like protein precursor (LOC100037730) carries an SNP that we consider the causal variegated phenotype mutation. Phylogenetic analysis shows the candidate is evolutionary closest to the Arabidopsis VAR1. The synonymous mutation created by the SNP generated a miRNA binding site, potentially disrupting the photoprotection mechanism and thylakoid development, resulting in leaf variegation. CONCLUSION: We described the histology, anatomy, physiology, and image analysis of four classes of cell layers and chloroplast degradation in a tomato plant with a variegated phenotype. The genomics and bioinformatics pipeline revealed a VAR1-related FtsH mutant, the first of its kind in tomato variegation phenotypes. The miRNA binding site of the mutated SNP opens the way to future studies on its epigenetic mechanism underlying the variegation.

Comparative analysis of chloroplast genomes of Pulsatilla species reveals evolutionary and taxonomic status of newly discovered endangered species Pulsatilla saxatilis.

BACKGROUND: Pulsatilla saxatilis, a new species of the genus Pulsatilla has been discovered. The morphological information of this species has been well described, but its chloroplast genome characteristics and comparison with species of the same genus remain to be reported. RESULTS: Our results showed that the total length of chloroplast (cp.) genome of P. saxatilis is 162,659 bp, with a GC content of 37.5%. The cp. genome contains 134 genes, including 90 known protein-coding genes, 36 tRNA genes, and 8 rRNA genes. P. saxatilis demonstrated similar characteristics to other species of genus Pulsatilla. Herein, we compared cp. genomes of 10 species, including P. saxatilis, and found that the cp. genomes of the genus Pulsatilla are extremely similar, with a length of 162,322-163,851 bp. Furthermore, The SSRs of Pulsatilla ranged from 10 to 22 bp in length. Among the four structural regions of the cp. genome, most long repeats and SSRs were detected in the LSC region, followed by that in the SSC region, and least in IRA/ IRB regions. The most common types of long repeats were forward and palindromic repeats, followed by reverse repeats, and only a few complementary repeats were found in 10 cp. genomes. We also analyzed nucleotide diversity and identified ccsA_ndhD, rps16_trnK-UUU, ccsA, and rbcL, which could be used as potential molecular markers for identification of Pulsatilla species. The results of the phylogenetic tree constructed by connecting the sequences of high variation regions were consistent with those of the cp. gene phylogenetic tree, and the species more closely related to P. saxatilis was identified as the P. campanella. CONCLUSION: It was determined that the closest species to P. saxatilis is P. campanella, which is the same as the conclusion based on pollen grain characteristics, but different from the P. chinensis determined based on morphological characteristics. By revealing information on the chloroplast characteristics, development, and evolution of the cp. genome and the potential molecular markers, this study provides effective molecular data regarding the evolution, genetic diversity, and species identification of the genus Pulsatilla.

Sequence characteristics, genetic diversity and phylogenetic analysis of the Cucurbita ficifolia (Cucurbitaceae) chloroplasts genome.

BACKGROUND: Curcubita ficifolia Bouché (Cucurbitaceae) has high value as a food crop and medicinal plant, and also has horticultural value as rootstock for other melon species. China is home to many different cultivars, but the genetic diversity of these resources and the evolutionary relationships among them, as well as the differences between C. ficifolia and other Cucurbita species, remain unclear. RESULTS: We investigated the chloroplast (cp) genomes of 160 C. ficifolia individuals from 31 populations in Yunnan, a major C. ficifolia production area in China. We found that the cp genome of C. ficifolia is ~151 kb and contains 128 genes, of which 86 are protein coding genes, 34 encode tRNA, and eight encode rRNAs. We also identified 64 SSRs, mainly AT repeats. The cp genome was found to contain a total of 204 SNP and 57 indels, and a total of 21 haplotypes were found in the 160 study individuals. The reverse repeat (IR) region of C. ficifolia contained a few differences compared with this region in the six other Cucurbita species. Sequence difference analysis demonstrated that most of the variable regions were concentrated in the single copy (SC) region. Moreover, the sequences of the coding regions were found to be more similar among species than those of the non-coding regions. The phylogenies reconstructed from the cp genomes of 61 representative species of Cucurbitaceae reflected the currently accepted classification, in which C. ficifolia is sister to the other Cucurbita species, however, different interspecific relationships were found between Cucurbita species. CONCLUSIONS: These results will be valuable in the classification of C. ficifolia genetic resources and will contribute to our understanding of evolutionary relationships within the genus Cucurbita.

Extraction and analysis of high-quality chloroplast DNA with reduced nuclear DNA for medicinal plants.

BACKGROUND: Obtaining high-quality chloroplast genome sequences requires chloroplast DNA (cpDNA) samples that meet the sequencing requirements. The quality of extracted cpDNA directly impacts the efficiency and accuracy of sequencing analysis. Currently, there are no reported methods for extracting cpDNA from Erigeron breviscapus. Therefore, we developed a suitable method for extracting cpDNA from E. breviscapus and further verified its applicability to other medicinal plants. RESULTS: We conducted a comparative analysis of chloroplast isolation and cpDNA extraction using modified high-salt low-pH method, the high-salt method, and the NaOH low-salt method, respectively. Subsequently, the number of cpDNA copies relative to the nuclear DNA (nDNA ) was quantified via qPCR. As anticipated, chloroplasts isolated from E. breviscapus using the modified high-salt low-pH method exhibited intact structures with minimal cell debris. Moreover, the concentration, purity, and quality of E. breviscapus cpDNA extracted through this method surpassed those obtained from the other two methods. Furthermore, qPCR analysis confirmed that the modified high-salt low-pH method effectively minimized nDNA contamination in the extracted cpDNA. We then applied the developed modified high-salt low-pH method to other medicinal plant species, including Mentha haplocalyx, Taraxacum mongolicum, and Portulaca oleracea. The resultant effect on chloroplast isolation and cpDNA extraction further validated the generalizability and efficacy of this method across different plant species. CONCLUSIONS: The modified high-salt low-pH method represents a reliable approach for obtaining high-quality cpDNA from E. breviscapus. Its universal applicability establishes a solid foundation for chloroplast genome sequencing and analysis of this species. Moreover, it serves as a benchmark for developing similar methods to extract chloroplast genomes from other medicinal plants.

Adaptive evolution of chloroplast division mechanisms during plant terrestrialization.

Despite extensive research, the origin and evolution of the chloroplast division machinery remain unclear. Here, we employ recently sequenced genomes and transcriptomes of Archaeplastida clades to identify the core components of chloroplast division and reconstruct their evolutionary histories, respectively. Our findings show that complete division ring structures emerged in Charophytes. We find that Glaucophytes experienced strong selection pressure, generating diverse variants adapted to the changing terrestrial environments. By integrating the functions of chloroplast division genes (CDGs) annotated in a workflow developed using large-scale multi-omics data, we further show that dispersed duplications acquire more species-specific functions under stronger selection pressures. Notably, PARC6, a dispersed duplicate CDG, regulates leaf color and plant growth in Solanum lycopersicum, demonstrating neofunctionalization. Our findings provide an integrated perspective on the functional evolution of chloroplast division machinery and highlight the potential of dispersed duplicate genes as the primary source of adaptive evolution of chloroplast division.

The dicot homolog of maize PPR103 carries a C-terminal DYW domain and may have a role in C-to-U editing of some chloroplast RNA transcripts.

In plants, cytidine-to-uridine (C-to-U) editing is a crucial step in processing mitochondria- and chloroplast-encoded transcripts. This editing requires nuclear-encoded proteins including members of the pentatricopeptide (PPR) family, especially PLS-type proteins carrying the DYW domain. IPI1/emb175/PPR103 is a nuclear gene encoding a PLS-type PPR protein essential for survival in Arabidopsis thaliana and maize. Arabidopsis IPI1 was identified as likely interacting with ISE2, a chloroplast-localized RNA helicase associated with C-to-U RNA editing in Arabidopsis and maize. Notably, while the Arabidopsis and Nicotiana IPI1 orthologs possess complete DYW motifs at their C-termini, the maize homolog, ZmPPR103, lacks this triplet of residues which are essential for editing. In this study we examined the function of IPI1 in chloroplast RNA processing in N. benthamiana to gain insight into the importance of the DYW domain to the function of the EMB175/PPR103/ IPI1 proteins. Structural predictions suggest that evolutionary loss of residues identified as critical for catalyzing C-to-U editing in other members of this class of proteins, were likely to lead to reduced or absent editing activity in the Nicotiana and Arabidopsis IPI1 orthologs. Virus-induced gene silencing of NbIPI1 led to defects in chloroplast ribosomal RNA processing and changes to stability of rpl16 transcripts, revealing conserved function with its maize ortholog. NbIPI1-silenced plants also had defective C-to-U RNA editing in several chloroplast transcripts, a contrast from the finding that maize PPR103 had no role in editing. The results indicate that in addition to its role in transcript stability, NbIPI1 may contribute to C-to-U editing in N. benthamiana chloroplasts.

From cyanobacteria and cyanophages to chloroplasts: the fate of the genomes of oxyphototrophs and the genes encoding photosystem II proteins.

Chloroplasts are the result of endosymbiosis of cyanobacterial organisms with proto-eukaryotes. The psbA, psbD and psbO genes are present in all oxyphototrophs and encode the D1/D2 proteins of photosystem II (PSII) and PsbO, respectively. PsbO is a peripheral protein that stabilizes the O2-evolving complex in PSII. Of these genes, psbA and psbD remained in the chloroplastic genome, while psbO was transferred to the nucleus. The genomes of selected cyanobacteria, chloroplasts and cyanophages carrying psbA and psbD, respectively, were analysed. The highest density of genes and coding sequences (CDSs) was estimated for the genomes of cyanophages, cyanobacteria and chloroplasts. The synonymous mutation rate (rS) of psbA and psbD in chloroplasts remained almost unchanged and is lower than that of psbO. The results indicate that the decreasing genome size in chloroplasts is more similar to the genome reduction observed in contemporary endosymbiotic organisms than in streamlined genomes of free-living cyanobacteria. The rS of atpA, which encodes the α-subunit of ATP synthase in chloroplasts, suggests that psbA and psbD, and to a lesser extent psbO, are ancient and conservative and arose early in the evolution of oxygenic photosynthesis. The role of cyanophages in the evolution of oxyphototrophs and chloroplastic genomes is discussed.

The chloroplast and mitochondrial genomes of two Gomphonema parvulum (Bacillariophyta) environmental isolates from South Carolina (United States) and Virginia (United States).

Gomphonema parvulum is a cosmopolitan freshwater diatom that is used as an indicator in water quality biomonitoring. In this study, we report the culturing of two geographically separated isolates from southeastern North America, their morphology, and the sequencing and assembly of their mitochondrial and chloroplast genomes. Morphologically, both strains fit G. parvulum sensu lato, but the frustules from a protected habitat in South Carolina were smaller than those cited in the historic data of this species from the same location as well as a second culture from Virginia. Phylogenetic analyses using the rbcL gene placed both within a clade with G. parvulum. Genetic markers, including full chloroplast and mitochondrial genomes and the nuclear small subunit rRNA gene region were assembled from each isolate. The organellar genomes of the two strains varied slightly in size due to small differences in intergenic regions with chloroplast genomes of 121,035 bp and 121,482 bp and mitochondrial genomes of 34,639 bp and 34,654 bp. The intraspecific pairwise identities of the chloroplast and mitochondrial genomes of these two isolates were 97.9% and 95.4%, respectively. Multigene phylogenetic analysis demonstrated a close relationship between G. parvulum, Gomphoneis minuta, and Didymosphenia geminata.

The newly assembled chloroplast genome of Aeluropus littoralis: molecular feature characterization and phylogenetic analysis with related species.

Aeluropus littoralis, a halophyte grass, is widely distributed from the Mediterranean to the Indian subcontinent through the Mongolian Gobi. This model halophyte has garnered increasing attention owing to its use as forage and its high tolerance to environmental stressors. The chloroplast genomes of many plants have been extensively examined for molecular, phylogenetic and transplastomic applications. However, no published research on the A. littoralis chloroplast (cp) genome was discovered. Here, the entire chloroplast genome of A. littoralis was assembled implementing accurate long-read sequences. The entire chloroplast genome, with an estimated length of 135,532 bp (GC content: 38.2%), has a quadripartite architecture and includes a pair of inverted repeat (IR) regions, IRa and IRb (21,012 bp each), separated by a large and a small single-copy regions (80,823 and 12,685 bp, respectively). The features of A. littoralis consist of 133 genes that synthesize 87 peptides, 38 transfer RNAs, and 8 ribosomal RNAs. Of these genes, 86 were unique, whereas 19 were duplicated in IR regions. Additionally, a total of forty-six simple sequence repeats, categorized into 32-mono, four-di, two-tri, and eight-tetranucleotides, were discovered. Furthermore, ten sets of repeats greater than 20 bp were located primarily in the LSC region. Evolutionary analysis based on chloroplast sequence data revealed that A. littoralis with A. lagopoides and A. sinensis belong to the Aeluropodinae subtribe, which is a sister to the Eleusininae in the tribe Cynodonteae and the subfamily Chloridoideae. This subfamily belongs to the PACMAD clade, which contains the majority of the C4 photosynthetic plants in the Poaceae. The newly constructed A. littoralis cp genome offers valuable knowledge for DNA barcoding, phylogenetic, transplastomic research, and other biological studies.

How Did Thylakoids Emerge in Cyanobacteria, and How Were the Primary Chloroplast and Chromatophore Acquired?

The emergence of thylakoid membranes in cyanobacteria is a key event in the evolution of all oxygenic photosynthetic cells, from prokaryotes to eukaryotes. Recent analyses show that they could originate from a unique lipid phase transition rather than from a supposed vesicular budding mechanism. Emergence of thylakoids coincided with the great oxygenation event, more than two billion years ago. The acquisition of semi-autonomous organelles, such as the mitochondrion, the chloroplast, and, more recently, the chromatophore, is a critical step in the evolution of eukaryotes. They resulted from primary endosymbiotic events that seem to share general features, i.e., an acquisition of a bacterium/cyanobacteria likely via a phagocytic membrane, a genome reduction coinciding with an escape of genes from the organelle to the nucleus, and, finally, the appearance of an active system translocating nuclear-encoded proteins back to the organelles. An intense mobilization of foreign genes of bacterial origin, via horizontal gene transfers, plays a critical role. Some third partners, like Chlamydia, might have facilitated the transition from cyanobacteria to the early chloroplast. This chapter further details our current understanding of primary endosymbiosis, focusing on primary chloroplasts, thought to have appeared over a billion years ago, and the chromatophore, which appeared around a hundred years ago.

Diversification of Plastid Structure and Function in Land Plants.

Plastids represent a largely diverse group of organelles in plant and algal cells that have several common features but also a broad spectrum of morphological, ultrastructural, biochemical, and physiological differences. Plastids and their structural and metabolic diversity significantly contribute to the functionality and developmental flexibility of the plant body throughout its lifetime. In addition to the multiple roles of given plastid types, this diversity is accomplished in some cases by interconversions between different plastids as a consequence of developmental and environmental signals that regulate plastid differentiation and specialization. In addition to basic plastid structural features, the most important plastid types, the newly characterized peculiar plastids, and future perspectives in plastid biology are also provided in this chapter.

The Use of Nanopore Sequencing to Analyze the Chloroplast Transcriptome Part II: Bioinformatic Analyzes and Virtual RNA Blots.

Nanopore sequencing of full-length cDNAs offers unprecedented details of the plastid RNA metabolism. After the generation of the nanopore reads, several bioinformatic steps are required to analyze the data. In this chapter, we describe in a few simple command lines the processing and mapping of the reads as well as the generation of virtual Northern blots as a simple and familiar way to visualize Nanopore sequencing data.

The Use of Nanopore Sequencing to Analyze the Chloroplast Transcriptome Part I: Library Preparation.

Global understanding of plastid gene expression has always been impaired by its complexity. RNA splicing, editing, and intercistronic processing create multiple transcripts isoforms that can hardly be resolved using traditional molecular biology techniques. During the last decade, the wide adoption of RNA-seq-based techniques has, however, allowed an unprecedented understanding of all the different steps of chloroplast gene expression, from transcription to translation. Current strategies are nonetheless unable to identify and quantify full length transcripts isoforms, a limitation that can now be overcome using Nanopore Sequencing. We here provide a complete protocol to produce, from total leaf RNA, cDNA libraries ready for Nanopore sequencing. While most Nanopore protocols take advantage of the mRNA polyA tail we here first ligate an RNA adapter to the 3' ends of the RNAs and use it to initiate the template switching reverse transcription. The cDNA is then prepared and indexed for use with the regular Oxford Nanopore v14 chemistry. This protocol is of particular interest to researchers willing to simultaneously study the multiple post-transcriptional processes prevalent in the chloroplast.

Structure of the multi-subunit chloroplast RNA polymerase.

Chloroplasts contain a dedicated genome that encodes subunits of the photosynthesis machinery. Transcription of photosynthesis genes is predominantly carried out by a plastid-encoded RNA polymerase (PEP), a nearly 1 MDa complex composed of core subunits with homology to eubacterial RNA polymerases (RNAPs) and at least 12 additional chloroplast-specific PEP-associated proteins (PAPs). However, the architecture of this complex and the functions of the PAPs remain unknown. Here, we report the cryo-EM structure of a 19-subunit PEP complex from Sinapis alba (white mustard). The structure reveals that the PEP core resembles prokaryotic and nuclear RNAPs but contains chloroplast-specific features that mediate interactions with the PAPs. The PAPs are unrelated to known transcription factors and arrange around the core in a unique fashion. Their structures suggest potential functions during transcription in the chemical environment of chloroplasts. These results reveal structural insights into chloroplast transcription and provide a framework for understanding photosynthesis gene expression.

Genome-wide identification of GA2ox genes family and analysis of PbrGA2ox1-mediated enhanced chlorophyll accumulation by promoting chloroplast development in pear.

BACKGROUND: Chlorophyll (Chl) is an agronomic trait associated with photosynthesis and yield. Gibberellin 2-oxidases (GA2oxs) have previously been shown to be involved in Chl accumulation. However, whether and how the PbrGA2ox proteins (PbrGA2oxs) mediate Chl accumulation in pear (Pyrus spp.) is scarce. RESULTS: Here, we aimed to elucidate the role of the pear GA2ox gene family in Chl accumulation and the related underlying mechanisms. We isolated 13 PbrGA2ox genes (PbrGA2oxs) from the pear database and identified PbrGA2ox1 as a potential regulator of Chl accumulation. We found that transiently overexpressing PbrGA2ox1 in chlorotic pear leaves led to Chl accumulation, and PbrGA2ox1 silencing in normal pear leaves led to Chl degradation, as evident by the regreening and chlorosis phenomenon, respectively. Meanwhile, PbrGA2ox1-overexpressing (OE) tobacco plants discernably exhibited Chl built-up, as evidenced by significantly higher Pn and Fv/Fm. In addition, RNA sequencing (RNA-seq), physiological and biochemical investigations revealed an increase in abscisic acid (ABA), methyl jasmonate (MeJA), and salicylic acid (SA) concentrations and signaling pathways; a marked elevation in reducing and soluble sugar contents; and a marginal decline in the starch and sucrose levels in OE plants. Interestingly, PbrGA2ox1 overexpression did not prominently affect Chl synthesis. However, it indeed facilitated chloroplast development by increasing chloroplast number per cell and compacting the thylakoid granum stacks. These findings might jointly contribute to Chl accumulation in OE plants. CONCLUSION: Overall, our results suggested that GA2oxs accelerate Chl accumulation by stimulating chloroplast development and proved the potential of PbrGA2ox1 as a candidate gene for genetically breeding biofortified pear plants with a higher yield.

Comparative analysis of chloroplast and mitochondrial genomes of sweet potato provides evidence of gene transfer.

The increasing number of plant mitochondrial DNA genomes (mtDNA) sequenced reveals the extent of transfer from both chloroplast DNA genomes (cpDNA) and nuclear DNA genomes (nDNA). This study created a library and assembled the chloroplast and mitochondrial genomes of the leafy sweet potato better to understand the extent of mitochondrial and chloroplast gene transfer. The full-length chloroplast genome of the leafy sweet potato (OM808940) is 161,387 bp, with 132 genes annotated, including 87 protein-coding genes, 8 rRNA genes, and 37 tRNA genes. The mitochondrial genome (OM808941) was 269,578 bp in length and contained 69 functional genes, including 39 protein-coding genes, 6 rRNA genes, and 24 tRNA genes. 68 SSR loci were found in the leafy sweet potato organelle genome, including 54 in the chloroplast genome and 14 in the mitochondria genome. In the sweet potato mitochondrial genome, most genes have RNA editing sites, and the conversion ratio from hydrophilic amino acids to hydrophobic amino acids is the highest, reaching 47.12%. Horizontal transfer occurs in the sweet potato organelle genome and nuclear genome. 40 mitochondrial genome segments share high homology with 14 chloroplast genome segments, 33 of which may be derived from chloroplast genome horizontal transfer. 171 mitochondrial genome sequences come from the horizontal transfer of nuclear genome. The phylogenetic analysis of organelle genes revealed that the leafy sweet potato was closely related to the tetraploid wild species Ipomoea tabascana and the wild diploid species Ipomoea trifida.

Comparative Analysis of Chloroplast Pan-Genomes and Transcriptomics Reveals Cold Adaptation in Medicago sativa.

Alfalfa (Medicago sativa) is a perennial forage legume that is widely distributed all over the world; therefore, it has an extremely complex genetic background. Though population structure and phylogenetic studies have been conducted on a large group of alfalfa nuclear genomes, information about the chloroplast genomes is still lacking. Chloroplast genomes are generally considered to be conservative and play an important role in population diversity analysis and species adaptation in plants. Here, 231 complete alfalfa chloroplast genomes were successfully assembled from 359 alfalfa resequencing data, on the basis of which the alfalfa chloroplast pan-genome was constructed. We investigated the genetic variations of the alfalfa chloroplast genome through comparative genomic, genetic diversity, phylogenetic, population genetic structure, and haplotype analysis. Meanwhile, the expression of alfalfa chloroplast genes under cold stress was explored through transcriptome analysis. As a result, chloroplast genomes of 231 alfalfa lack an IR region, and the size of the chloroplast genome ranges from 125,192 bp to 126,105 bp. Using population structure, haplotypes, and construction of a phylogenetic tree, it was found that alfalfa populations could be divided into four groups, and multiple highly variable regions were found in the alfalfa chloroplast genome. Transcriptome analysis showed that tRNA genes were significantly up-regulated in the cold-sensitive varieties, while rps7, rpl32, and ndhB were down-regulated, and the editing efficiency of ycf1, ycf2, and ndhF was decreased in the cold-tolerant varieties, which may be due to the fact that chloroplasts store nutrients through photosynthesis to resist cold. The huge number of genetic variants in this study provide powerful resources for molecular markers.

The Complete Mitochondrial Genome of Paeonia lactiflora Pall. (Saxifragales: Paeoniaceae): Evidence of Gene Transfer from Chloroplast to Mitochondrial Genome.

Paeonia lactiflora (P. lactiflora), a perennial plant renowned for its medicinal roots, provides a unique case for studying the phylogenetic relationships of species based on organelle genomes, as well as the transference of DNA across organelle genomes. In order to investigate this matter, we sequenced and characterized the mitochondrial genome (mitogenome) of P. lactiflora. Similar to the chloroplast genome (cpgenome), the mitogenome of P. lactiflora extends across 181,688 base pairs (bp). Its unique quadripartite structure results from a pair of extensive inverted repeats, each measuring 25,680 bp in length. The annotated mitogenome includes 27 protein-coding genes, 37 tRNAs, 8 rRNAs, and two pseudogenes (rpl5, rpl16). Phylogenetic analysis was performed to identify phylogenetic trees consistent with Paeonia species phylogeny in the APG Ⅳ system. Moreover, a total of 12 MTPT events were identified and 32 RNA editing sites were detected during mitogenome analysis of P. lactiflora. Our research successfully compiled and annotated the mitogenome of P. lactiflora. The study provides valuable insights regarding the taxonomic classification and molecular evolution within the Paeoniaceae family.

CLAW: An automated Snakemake workflow for the assembly of chloroplast genomes from long-read data.

Chloroplasts are photosynthetic organelles in algal and plant cells that contain their own genome. Chloroplast genomes are commonly used in evolutionary studies and taxonomic identification and are increasingly becoming a target for crop improvement studies. As DNA sequencing becomes more affordable, researchers are collecting vast swathes of high-quality whole-genome sequence data from laboratory and field settings alike. Whole tissue read libraries sequenced with the primary goal of understanding the nuclear genome will inadvertently contain many reads derived from the chloroplast genome. These whole-genome, whole-tissue read libraries can additionally be used to assemble chloroplast genomes with little to no extra cost. While several tools exist that make use of short-read second generation and third-generation long-read sequencing data for chloroplast genome assembly, these tools may have complex installation steps, inadequate error reporting, poor expandability, and/or lack scalability. Here, we present CLAW (Chloroplast Long-read Assembly Workflow), an easy to install, customise, and use Snakemake tool to assemble chloroplast genomes from chloroplast long-reads found in whole-genome read libraries (https://github.com/aaronphillips7493/CLAW). Using 19 publicly available reference chloroplast genome assemblies and long-read libraries from algal, monocot and eudicot species, we show that CLAW can rapidly produce chloroplast genome assemblies with high similarity to the reference assemblies. CLAW was designed such that users have complete control over parameterisation, allowing individuals to optimise CLAW to their specific use cases. We expect that CLAW will provide researchers (with varying levels of bioinformatics expertise) with an additional resource useful for contributing to the growing number of publicly available chloroplast genome assemblies.